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The 2D "Swiss-cheese" model of conducting media with round insulator inclusions is 
studied in the 2nd order of inclusion concentration and near the percolation threshold. 
The electric field distribution function is found to have power asymptotics for fields much 
exceeding the average field, independently on the vicinity to the threshold, due to finite 
probability of arbitrary proximity of inclusions. The strong field in the narrow necks 
between inclusions results in the induced persistent anisotropy of the system. The critical 
index for noise density is found, determined by the asymptotics of electric field distribution 
function. 



Introduction 

The distribution function of electric field is an important characteristics of a con- 
ducting mixture. It determines such properties as noise, electric breakdown and 
nonlinear dielectric permeabilities. Its long-range tails are usually considered as a 
result of the electric field amplification due to the vicinity to the percolation thresh- 
old i, 0,0. 

The purpose of the present study is to find the distribution function for the elec- 
tric field F in the conducting phase. The 2D model of conducting media containing 
randomly distributed round voids with density n and radius a is considered. This 
system is interesting both as a real object of nanoelectronics M and as an exact 
solvable model. Without inclusions the electric field is uniform. The field around a 
solitary smooth void is limited, so it does not produce the field, much exceeding the 
average field. The strong field appears in the 2nd order of void concentration only. 

We shall study both 2nd order of concentration and the field distribution at the 
percolation threshold. 

Potential distribution around two holes. 

The distribution of field in a conducting media is described by the solution of the 
continuity equation Vj = for current j = ctF = aVy?, where a is the conductivity 
of the conductor. The boundary conditions on the void boundaries require vanishing 
of the normal component of current density crnVy? (n is the normal to the circle). 
The electrostatic problem of field distribution is two-dimensional. 

Since the asymptotics of the field distribution is determined by narrow necks 
between two holes, we turn our attention on the media , containing two round holes 
only (Fig. 1). We assume that the x axis goes through circles centers. The conformal 
mapping of system of Fig.l to the band 1/Rq < Rew < Ro is given by the function 
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W = —b> (1) 
where z = x + iy, 25 is the distance between the holes with radii a, 

y/2a + 5 + V5 



b = V2a5 + S 2 , R 
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We should solve the problem of current flow with conditions ip ~ Im (F z) for 
z — > oo, where Fq is the average electric field. The corresponding potential in 
w-plane is the dipole potential with asymptotic behavior ip ~ Im (^f), for w — > 0. 

The dipole potential in the ring between two circles is determined by the sum of 
its images in circles: 
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The strong field results from the close approaching of two circles, i < a. In this 
limit (2) transforms: 

ip = nF a Im(cot — ) (3) 

z 

The last extreme limit means the absence of percolation between contacting circles 
if 5 — > 0. This leads to the gap of potential between up and down sides of the circles 
near the points of tangency: 

U = 2irF a (4) 

If the width of constriction is small but finite, the potential gap is distributed 
along the length of constriction \/2aS. The field in the vicinity of constriction is 

F = V-M 1 (5) 
y TT V a2S + y ya { ) 



The maximal value of field F max (5) = F ^2a/5 diverges with 5 — * 0. Another 
component of field, F x , remains finite if 5 — > and is insufficient below. 



Breakdown-induced training. 

The divergency of field and current density can lead to overheating and melting of 
narrow necks. If we assume that breakdown is determined by the critical density of 
heat emission Q c , all places where aF 2 > Q c will be overheated. It means that some 
number of overheated places occurs in any weak field F independently on the field 
strength. So the application of field irreversibly change the electrical properties of 
system. Moreover, the direction of field tends to the residual anisotropy of a system. 

We shall consider this effect in the framework of swiss-cheese system with rare, 
randomly distributed round voids. The residual change of linear conductivity in 
a low field originates from the breakdown of narrowest necks. The criterium of 
breakdown is F max (5) > F c = JQJa. 
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The contribution of necks between pairs of voids which have x-orientation and 
widths between 5 and S + dS to the average conductivity is 




Aa yy = 16nV2(na 2 ) 2 a J , (6) 



One should average (6) over 5 and the direction of pair of voids. The value 6 is limited 
by the maximum S max determined by the situation when the field in the neck exceeds 



the breakdown field F c = F y - n ^ !L cos9, where 9 is the angle between the external 
field and the direction of axis, connecting centers of circles. The resulting change of 
conductivity tensor is 

128>/2, 2n2/I F o|\ 3 



Aa xx = -^-(^T(^f) (7) 
Aa yy = AAa xx . (8) 

The direction y is determined along the training field. 

The change of longitudinal conductivity appears to be 4 times higher than 
transversal. The equation (8) shows the non-analytical dependence on the field. 
It should be noted that the training, being a nonlinear effect in applied field, never- 
theless means the change of linear conductivity. 

Distribution function of field for pairs of circles. 

The asymptotics of distribution function for high field is determined by rare ap- 
proaching of circles, when the distance 6 between them is much less than the radius 
and mean distance (7m) _1 / 2 . 

The probability of field to have value F, P$jj(F) for the fixed distance 5 and 
potential gap U is determined by area where F is between F and F + dF. By means 
of (H) we find: 

P W(F) = ^) 3/4 ^7i j p5/2 (f/ _ F7rv ^ ) l/2 

This result should be integrated by 5 with the density of pairs 4nn 2 a d8, distance 
of which is within the range 25 and 2(6 + d5) . The result of averaging on 5 and the 
direction of axis y is 

P(F) = ^(W) 2 |§ (10) 

The equation (10) gives power long-range tail of the electric field distribution 
function P(F) ~ F~ 6 . 

The long-range distribution tail gives the divergency of high order power of the 
local field, starting from 6th power and does not affect the density of 1/ f noise, 
determined by 4th power of F. 

The existence of power tail results from the microgeometry and is not connected 
with the vicinity to the percolation threshold. Thus such power tails are universal 
property of 2D system both far and near threshold. They result from the possibility 
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for 2D system to have a finite potential gap U on infinitely small distance, connected 
with the necessity of current to flow around tangent circles on the finite distance 
to equilibrate the potential gap. On the other hand, the shortest path between two 
infinitely close points along the conducting phase in the 3D case is always small and 
the field near close pair of inclusions is limited. Hence, no long-range tail exists in 
3D system of non-intersecting spheres and only exponential tail P(F) is possible in 
the intersecting inclusions 2D or 3D models determined by rare large coupling of 
inclusions. 



Percolational situation. 

The power-law tails due to formula (|10D are determined by very small part of total 
current while most part of current flows around both disks. Hence in the perco- 
lational limit na 2 ~ 1 the contribution retains, described by random necks inside 



the conducting phase. Formula (fLOf) may be applied to this case also (except for 
numerical coefficient) if one changes na 2 — ► 1. 

The key bonds tend to additional possibility for local field amplification near the 
threshold. 

The distribution function of voltages in lattice models was studied on the metal 
side of percolation threshold 0, §. It was found that this function has a long- 
range tail, P(U) ~ exp((ln(U / Uq)) 2 / C) , which is much weaker than any power-like 
one. The distribution function in insulator phase, found below, exhibits power 
behavior. In the Swiss-cheese model narrow necks, insufficient for current flow may 
be considered as almost insulating phase, so these tails are reproducing themselves 
in the contribution from narrow necks to metal phase distribution function long tail. 

To estimate the insulator phase distribution function tail we shall use the evident 
formula for energy: 

87T 87T 

where e is the permittivity of the insulator. Just below percolation threshold the 
effective permittivity e e diverges like r _<? , producing divergency of (F 2 ) if the system 
comes to the percolation threshold r = 0, while the mean field (F) is limited. These 
two facts do not contradict if the distribution function behaves like 1/F' 3 ~ r , where 
< r < 1. 

Near the percolation threshold one should expect the scaling behavior of a dis- 
tribution function, both above and below the threshold. The characteristic scale for 
voltage on single bond U = Fa is determined by the maximal voltage U m in a per- 
colation cell. This voltage can be estimated by assumption, that full voltage across 
the cell with size L c is applied to one disconnected bond: U m ~ F L C ~ F aT^ u . 
Hence the distribution function for voltages is 

P(U) = -^(F a/U) 3 - r f(^rn (12) 
rod r a 



The expression for exponent r = q/v is determined by equation (|iT|). According 
to known values of critical exponents in 2D case r = 0.9. This result is consistent 
with numerical result of [0] but slightly deviates from estimation based on the as- 
sumption, that asymptotics p(U) is determined by the number of single-disconnected 
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bond per a percolation cell || N sd ~ r -1 : p(U m ) ~ uh ~ T?>u ~ l - Like so called 
nodes-links-blobs model we shall disregard small deviations from 1 of critical expo- 
nents for effective conductivity a e ~ r* and effective permittivity £ e ~ t~ 9 . 

The far tail for electric field distribution function in the Swiss-cheese model on 
the metal side of percolation transition is determined by narrow necks on the infinite 
cluster. The strongest voltage is expected from a narrow neck on a branch of infi- 
nite cluster clumped to the backbone in two points with distance L c between them. 
The width of the neck is supposed so small, that it resistance exceeds the backbone 
resistance and current through it is negligible. Hence the typical width is limited 
by 5 < r 2t . Such necks may be considered as insulating for the problem of current 
flow. For the neck with the potential gap has order of magnitude U = FqL c . Inte- 
grating the distribution function of field ([D together with the distribution function 
of voltages we find the distribution function of field in the vicinity of percolation 
threshold 



P(F) ~r-('+ 3 ^^ for F>F T- {t+u) (13) 
F 

^3/2 

P(F) ~ r^' 2 ^ for F < F r^ t+V \ (14) 

Comparison of formulae (0), ([14]) and ( |T0D shows that the validity of (|T4D is lim- 
ited by an inequality FoT~( 5t+u ^ 7 < F < Fqt~^ +v ' where the contribution of single- 
disconnected bonds prevails. For lower fields F < F r^^ 5t+U ^ 7 the asymptotics of 
distribution function is given by the contribution of narrow necks inside metal phase 

(0)- 

The amplification of field distribution function tails near the percolation thresh- 
old leads to enhancement of current noise. The integral noise in disordered media 
is determined by mean 4th power of field, according to |J: 

c - mwtf- ( ' 



The substitution of (|T3| , |14[) to (15) gives the critical exponents for integral noise: 



Concluding remarks. 

Let us summarize the obtained results. 

We found the field distribution around pair of circular insulating inclusions into 
the conducting phase. It was used to find the second virial correction on the density 
of inclusions to the distribution function of field. The distribution function was found 
to have the power-like long range tail, determined by close locations of inclusions. 
The same configurations result in the training of medium, that is, an irreversible 
modification of material by a weak electric field, which destroys the narrow necks 
and produces the persistent anisotropy. 

It was demonstrated, that the tail of distribution function is enhanced near 
the percolation threshold. The field distribution function was employed to find 



5 



the current noise in the system. It was found that it is amplified of the near the 
percolaton threshold. The main sources of this amplification are the narrow necks 
belonging to the percolation cluster backbone. 

The results of the present article have close tights to the system dimensional- 
ity. They are determined by the possibility of existence of a finite potential gap 
between two close points in 2D, in contradiction with 3D system. This conclusion 
is independent on the model of inclusions we considered. 

The work was partially supported by Russian Foundation for Basic Researches 
(Grants 950204432 and 960219353) and Folkswagen-Stiftung. 
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